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Abstract 


Potential curves for the lowest two electronic states (X 2 E^ and A 2 n„) of 
Al^ have been computed using complete active space SCF/multireference Cl wave 
functions and large Gaussian basis sets. The lowest observable vertical ionization 
potential (to Alj -X' 2 E^ ) of the Al 2 A 3 II U ground state is calculated to occur 
around 6.1 eV, in excellent agreement with the experimental range of 6.0-6.42 eV 
obtained in recent cluster ionization studies by Cox and co-workers. The second 
vertical ionization potential (to Al^ A 2 II U ) occurs near 6.4 eV, also within the 
experimental range. The adiabatic IP of 5.90 eV is in good agreement with the value 
of 5.8-6. 1 eV deduced by Hanley and co-worker: from the difference in thresholds 
between collision induced dissociation processes of AI3 . The computed IP values are 
somewhat larger than those deduced from bran cling ratios in cluster fragmentation 
experiments by Jarrold and co-workers. The observation of an ionization threshold 
below 6.42 eV is shown to be incompatible wr h an Al 2 ground electronic state 
assignment of 3 E“, but the separation between the two lowest states of Al 2 is so 
small that it is likely that both are populated in the experiments, so that this 
does not provide unambiguous support for the recent theoretical assignment of the 
ground state as 3 II U . 


(AASA-IH-101171) 
iClEBUAlS Cf A 12 


I BE LOW ESI 2CBI2A1LOH 

(NASA) 1b F ^ scl 20H 


G 3/72 


N69-106C5 


Unclas 

0158953 


f Mailing address: NASA Ames Research Center, Moffett Field, CA 94035 


1 



I. Introduction 

Recent studies of aluminum clusters by J< rrold et al. [l], Hanley et al. [2], 
and Cox et al. [3] have provided estimates of cluster ionization potentials. While 
these sources agree for most ionization potenti ils of Al n , n < 15, Jarrold et al. 
inferred a value for AI 2 of 5.5 eV, outside the range of 6.0-6.4 eV obtained by Cox 
et al. Hanley et al. deduced a value of 5.8 eV, but noted that the dissociation 
might involve a spin-forbidden process; correcting for this gave a value of 6.2 eV. 
(This should be revised to 6.1 eV using the splitting between the lowest singlet and 
triplet states of AI 2 (a 1 !!^ - A 3 n u ) from the best available calculation [4]). The 
work of Cox et al. involves laser photoionization of a cluster beam followed by mass 
spectrometric analysis of the resulting cluster ions, and it would be expected that 
the failure to observe production of Al^" with laser energies below 6 eV provides a 
reliable lower bound for the ionization potential of AI 2 . The result of Jarrold et al. 
is obtained from analysis of the branching ratios in observations of fragmentation of 
Al^ cluster ions, and is therefore a less direct measurement than that of Cox et al. 
Using the data analysis method of Jarrold et a/., Hanley et al. obtained a value of 
5.2 eV. Upton [5] has calculated a vertical IP for the -Y 3 H U state of AI 2 of 6.02 eV. 
However, the same level of theory incorrectly predicts the ground state of AI 2 to be 
3 E“, lying 0.06 eV below the 3 n^ state [6]. Sunil and Jordan [7] have computed 
vertical and adiabatic IPs of 6.12 and 5.92 eV using the coupled-cluster double 
excitation approach including single and triple excitations (CCD-)-ST(CCD)). This 
approach was in reasonable agreement with the t est calculations for the AI 2 A 3 n u - 
A 3 £~ separation [4] and the calculated IP is in good agreement with the work of 
Cox et al. However, since the available experimental values for the AI 2 IP vary so 
markedly, it is worthwhile reconsidering this pre blem with a higher level of theory 
than has been used previously. 

In a recent study of the low-lying states cf AI 2 , Bauschlicher et al. [4] con- 
cluded that the ground state was 3 n u , but that the first excited state, 3 E~, was 
only 180 cm -1 higher in energy, with the state lying some 2 000 cm -1 above 
the ground state. While the separation between the lowest two states is very small, 
the reliability of the prediction was established by performing full Cl benchmark 
calculations to calibrate the appropriate comp Rational method (a complete ac- 
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tive space SCF (CASSCF) calculation followed by a multireference Cl calculation 
with all CASSCF configurations used as reference configurations — second-order 
Cl (SOCI)); this method was then applied in a ' ery large atomic basis set, so that 
errors due to basis set incompleteness or an inadequate treatment of correlation 
were essentially eliminated. As this SOCI approach has been shown to be very 
reliable in the calculation of ionization potentials (IPs) of other systems with the 
same number of valence electrons as Al 2 [8], this approach offers a way of accurately 
computing the IPs of Al 2 , and thereby characterizing the low-lying states of Al^" 
and possibly rationalizing the discrepancies betveen the existing experiments. 

In the next section we describe briefly th< computational methods (a more 
complete discussion is given in Ref. 4), in sec ion III we discuss the computed 
results, and our conclusions are given in section IV. 

II. Methods 

The contracted Gaussian basis used is the [ is hp 3 d 2f] set from Ref. 4, which 
is a general contraction, based on atomic natural orbitals, of a (20.s 13p 6 d 4 /) prim- 
itive basis. All results presented in this work we re obtained using CASSCF/SOCI 
wave functions. The active space for the CASSCF calculations comprised those 
molecular orbitals derived from the atomic 3s and 3 p orbitals, with six active elec- 
trons for the neutral and five for the ionic states. Only these electrons were cor- 
related in the SOCI calculations. This procedure is known to produce excellent 
spectroscopic constants for Al 2 , where it has been calibrated against full Cl cal- 
culations, and it can be expected to perform equally well for Al^\ In this regard, 
we note that for r(Al-Al)=100 ao, the IP is 5.93 eV, or 0.05 eV smaller than the 
experimental value for A1 atom [9]. Thus we expect our computed IPs of Al 2 to be 
systematically about 0.05-0.10 eV too small. 

Spectroscopic constants (r e , u; e ) were obtained from a fit in 1/r to three points 
around the minimum of each curve. D e values w >re obtained using a supermolecule 
approach at r = 100 a 0 to reduce size-consister cy errors. For the calculation of 
Franck-Condon factors, the potential curves weie represented by a tensioned cubic 
spline interpolation of the computed energy points, and vibrational wave functions 
were obtained by numerical integration. 
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All calculations were performed using th< MOLECULE-SWEDEN [10,11] 
codes on the NAS Facility CRAY 2. 

III. Results and discussion 

Potential curves for the A 3 II U and A 3 ?,' stales of Al 2 and the A 2 E+ and A 2 H U 
states of Al^ are displayed in Fig. 1, and the corresponding spectroscopic constants 
are listed in Table I. The Al 2 A 3 II U ground state wavefunction is dominated by the 
open-shell configuration < 7 ^, while the A 3 Ej etcited state wavefunction is domi- 
nated by the open-shell configuration 7r 2 . Since f shorter bond distance is required 
for effective t bonding, the bond length of the A 3 E U state would be expected to 
be longer, as observed, and the same effect manifests itself in the Al 2 states: the 
A 2 £+ ground state has a a g bonding electron and a considerably longer bond than 
the A 2 n u state with its single 7r u bonding elec ron. Although the bond lengths 
in all the neutral and ionic states listed in Table I are rather large, all the states 
have significant (% 1 eV or more) binding energies. It is interesting to note that for 
Al 2 , the difference in energy between the A 3 E S and the A 3 II U states is very small 
(0.02 eV), whereas the difference in energy for the A 2 E^ and A 2 II U states of Al 2 is 
much larger (0.52 eV). This is despite the fact that the occupations of the Al 2 and 
AlJ states differ only in the single 7r u electron. The reason for the large difference 
comes from the increased electrostatic contribution to the bonding in Al 2 ■ The 
<r g electron in the 2 £+ state shields the positively charged Al centers from each 
other, thus reducing the electrostatic repulsion this effect is not present for the 
2 II U state. In addition, the much shorter bond length of the 2 II U state leads to a 
greatly increased electrostatic repulsion compart d to the 2 E^ state. Thus the 
state is significantly more stable than the 2 I1 U s ate for Al^. Finally, we note that 
the spectroscopic constants of Table I are in good agreement with those obtained 
by Sunil and Jordan [7] using a smaller basis set and the CCD+ST(CCD) method. 

It is clear from the potential curves of Fig. 1 that the A 3 !!*— * A 2 n u ionization 
will have very similar values for the adiabatic and vertical IPs, whereas the large 
difference in r e for the ground and ionized states :!or the A 3 !!*— ► A r2 E^ and A 3 H g — + 
A 2 n u ionizations leads to a considerable difference between their adiabatic and 
vertical IPs. The IP values listed in Table II confirm these qualitative deductions. 
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We note that since the A 3 £j state cannot ionize to the J¥ 2 E^ state by a one- 
electron process, the intensity for this ionizatioi is expected to be very low. Like 
the spectroscopic constants, the IPs are in good agreement with the results of Sunil 
and Jordan [7]. The vertical values correspond to that ionization from the v" = 0 
level of the appropriate neutral dimer state that has the largest Franck- Condon 
factor, as given in Table III. Only ionizations f om v" = 0 levels are listed: other 
values may be calculated from the potential cur ves, the computed energies defining 
which are listed in the Appendix. The most probable transitions for the lowest 
ionization process, .A 3 II U — > .A 2 £+, are actually to the levels around v' = 12 at 
49 537 cm -1 , some 2000 cm -1 above the adiab itic value. The next lowest energy 
ionization process is -A 3 II U — ► A 2 II U , which has a vertical IP of around 51 851 cm -1 , 
the same as the adiabatic value. Finally, ionization from the A 3 E“ excited state will 
occur mainly to the levels v' = 2 and somewhat higher, beginning at 52 080 cm -1 . 

These results are completely consistent with the experimental observations of 
Cox et al. [3] We predict that there w T ill be litile or no intensity from ionization 
into the lowest levels of Al^ , so that although the lowest adiabatic IP for 

AI 2 is below 6 eV (and is, in fact, some 600 cm -1 less than that of Al atom), the 
experimental threshold for observing ionization of AI 2 will indeed be between 6.0 
and 6.42 eV, as noted by Cox et al. At photo onization energies above 6.42 eV, 
ionization from the AI 2 .A 3 !!,, state to give the A 2 II U state of the ion will also be 
observed. Further, if (as seems very likely under the experimental conditions of 
Ref. 3) the A 3 E~ state of AI 2 is significantly populated, ionization from this state 
to the A 2 n u ion state will also occur at photon energies above 6.42 eV. 

In addition to providing a complete explanation of the ionization threshold 
behavior observed by Cox et al. [3], the above Jesuits suggest another method for 
confirming the assertion in Ref. 4 that the ground state of AI 2 is indeed 3 n v . While 
the full Cl calibration calculations described in Ref. 4 should provide a very small 
error bar for the prediction of a 3 II u ground-stale, it must be admitted that the es- 
timated separation between the two lowest triple: states is very small. However, the 
3 £j state can only ionize to give 2 II U , a process v ith a computed energy threshold of 
almost 6.5 eV, some 0.4 eV above the onset of ionization from 3 D^. The error asso- 
ciated with this theoretical prediction can hardly be greater than a few hundredths 
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of an eV. Therefore, if Al 2 molecules could be generated in a source that allowed 
almost complete cooling of the electronic degree of freedom, the observation of an 
ionization threshold of around 6 eV would provide absolute confirmation of a 3 II U 
ground state for Al 2 . However, we note that th< presence of A1 atoms, which also 
ionize above this threshold, may complicate performing such an experiment. 

While the current results agree well with Ihe experimental results of Cox et 
ah [3], they differ considerably from those of Ja rold et ah [ 1 ] These workers infer 
an (essentially adiabatic) ionization potential of i .round 5.5 eV for Al 2 , considerably 
smaller than our estimate. However, the analysi of cluster ion fragmentation used 
to obtain this estimate in Ref. 1 may not be r el able for Al^, which is the species 
used to estimate the Al 2 IP. As noted previously, Hanley et ah [2] obtain a value 
of only 5.2 eV using the analysis of Jarrold et al , compared with 5.8 eV from their 
own analysis. While cluster ions A1+ with n < 1 1 primarily fragment to Al n _i and 
A1+ , Al^ fragments to give mainly Al^" and Al. This in turn suggests that a number 
of the hypotheses used in inferring the ionizatio i potentials, such as the similarity 
between transition states on various dissociatio i pathways, may not be valid. It 
should be noted that Jarrold et al. also point to the fact that dissociation of Al^ 
to give Al 2 and Al + is spin-forbidden as a possible source of error in their analysis 
(or at least in comparison of their work with that of other groups). However, as 
noted above, this is not observed to be the primary fragmentation pathway for AI 3 , 
and the observed fragmentation to give Al^" anc Al is spin-allowed independent of 
whether the ground state of Al^~ is a triplet or a singlet. It seems improbable that 
this is the basis for the disagreement between the results of Jarrold et ah w r ith those 
of Cox et ah and the present work: it seems mere likely that the model of cluster 
ion fragmentation used in Ref. 1 is simply not v.Jid for the smallest cluster ions. 

The result of Hanley et ah [ 2 ] (5.8 eV) is expected to correspond to the adiabatic 
IP. As they noted, the Al^“ — ►Al 2 -fAl + path is expected to be spin forbidden and 
therefore their IP should be corrected by the splitting between the lowest singlet and 
triplet states for Al 2 . They used the separation of 0.4 eV from the work of 

Upton [ 6 ], which incorrectly predicts a 3 £” ground state for Al 2 . However, better 
calculations [4] yield a value of 0.3 eV for the singlet-triplet correction, based on 
the 3 £+- 3 n u separation; this would result in an IP of 6.1 eV. This is in reasonable 
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agreement with our computed value of 5.90 eV (~5. 95-6. 00 eV if we correct our 
IP for the error in the computed A1 atom IP). T rus it appears that the method of 
Hanley et al. obtains reasonable estimates for the cluster IPs. It should be noted 
that their method should become more reliable is the clusters increase in size, as 
the difference in energy between states of different spin will be even smaller. 

Our computed Al 2 adiabatic IP of 5.90 eV is slightly less than the Al atom 
IP of 5.93 eV [9], using the same level of treatment. Thus the Al 2 IP is smaller 
than the Al IP, as suggested by Jarrold et al. [1], although only by a very small 
amount. Upton [5] has computed a vertical IP for Al 2 of 6.02 eV for the 3 n u 
state, in a smaller basis set and at a lower leve of correlation treatment than in 
the present work. As we noted above this level of treatment incorrectly predicts a 
3 E - ground state for Al 2 , thus the modest diffe ence with the present work is not 
unexpected. The results of Sunil and Jordan [<1 are in good agreement with the 
present calculations, and therefore we conclude that the CCD-I-ST(CCD) level of 
treatment gives a good description of this system. 


IV. Conclusions. 

Large-scale multireference Cl calculations in an extended basis set predict the 
lowest adiabatic IP of Al 2 to be 5.90 eV, very close to that of Al atom. However, 
the lowest vertical IP is predicted to be at least 0.25 eV larger than the adiabatic 
value because of the large increase in r e on ioni 'ation. This is consistent with the 
observation of Cox et al. [3] that the threshold for ionization is between 6.0 and 
6.42 eV, and suggests that the inferred adiabat:c IP of 5.5 eV of Jarrold et al. [1] 
is much too low. The IP of Hanley et al. [2], if < orrected with the best estimate of 
the singlet-triplet separation in Al 2 , is also in g.x>d agreement with the computed 
values. At ionization energies above 6.4 eV, ionization to give the excited Al^ A 2 H^ 
state can occur, as can ionization from the A 3 1 ~ state of Al 2 if the experimental 
conditions are such that this state is appreciably populated. It is suggested that 
the observation of ionization below 6.42 eV in electronically cold Al 2 would confirm 
the assertion that the ground state of Al 2 is indeed 3 n„. 
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Table I. Al£ and Alj spectroscopic constants. 



r e (A) 

u? e (cm 1 ) 

D e (eV) 

T e (cm 

X 3 U u 

2.727 

278 

1.40 


a 3 E; 

2.492 

345 

1.38 

165 


3.212 

169 

1.42 


A 2 n„ 

2.794 

218 

0.90 

4214 


° From Ref. 4. 
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Table II. Adiabatic and vertical ionization potentials (eV) of AI 2 . 



Adiabatic 

Vertical 

x 3 u u - 

5.90 

6.14 

A 3 E“ - A 2 n u 

6.40 

6.53 

A 3 n u - A 2 n u 

6.43 

6.43 
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Table III. Franck-Condon factors (<7r'o) for ionizition of AI 2 from t>”=0. 



X 3 U u 

- a 2 e+ 


-4 A 2 ll u 

A 3 n„ 

- A 2 Jl u 

v' 

Qv’o 

A£(cm- 3 ) 

9v'o 

A£’(cm“ 1 ) 

9r'0 

AF(cm _1 ) 

0 

0.0000 

47 612 

0.7855 

5: 851 

0.0065 

51 651 

1 

0.0003 

47 781 

0.1933 

51’ 067 

0.0268 

51 868 

2 

0.0012 

47 948 

0.0204 

51’ 281 

0.0591 

52 081 

3 

0.0036 

48 114 

0.0009 

51’ 490 

0.0937 

52 290 

4 

0.0081 

48 278 



0.1196 

52 496 

5 

0.0152 

48 441 



0.1305 

52 698 

6 

0.0250 

48 602 



0.1264 

52 897 

7 

0.0367 

48 761 



0.1114 

53 092 

8 

0.0491 

48 919 



0.0912 

53 283 

9 

0.0608 

49 076 



0.0704 

53 472 

10 

0.0703 

49 231 



0.0519 

53 657 

11 

0.0769 

49 385 



0.0370 

53 838 

12 

0.0799 

49 537 



0.0256 

54 017 

13 

0.0795 

49 688 



0.0173 

54 192 

14 

0.0761 

49 837 



0.0115 

54 364 

15 

0.0704 

49 985 



0.0076 

54 532 

16 

0.0632 

50 132 



0.0049 

54 697 

17 

0.0552 

50 277 



0.0032 

54 859 

18 

0.0470 

50 421 



0.0020 

55 018 

19 

0.0392 

50 564 



0.0013 

55 173 

20 

0.0321 

50 705 



0.0008 

55 325 

21 

0.0258 

50 845 



0.0005 

55 473 

22 

0.0204 

50 984 



0.0003 

55 618 

23 

0.0159 

51 122 



0.0002 

55 759 

24 

0.0123 

51 258 



0.0001 

55 897 
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Appendix. The total SOCI energies in Eh, relative to -483.0. 
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-0.527 204 
-0.628 335 


-0.662 954 

-0.646 938 

-0.671 830 

-0.664 080 

-0.677 145 

-0.676 493 

-0.679 961 
-0.680 682 

-0.685 345 

-0.681 049 
-0.681 125 

-0.691 514 

-0.680 962 

-0.695 658 

-0.680 094 

-0.698 269 

-0.678 720 

-0.699 721 
-0.700 103 
-0.700 297 
-0.700 326 

-0.675 181 

-0.700 213 

-0.673 248 

-0.699 633 


-0.698 684 

-0.669 398 

-0.697 462 
-0.696 040 
-0.694 474 
-0.692 810 

-0.663 371 

-0.691 080 
-0.689 312 

-0.659 859 
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-0.900 386 

-0.885 410 

-0.892 405 

-0.878 458 

-0.885 566 

-0.873 790 

-3.880 113 

-0.870 842 


-0.866 026 
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Figure Caption. 

The AI2 and Alj potential energy curves, relat ve to the AI2 A^ 3 II U ground state 
minimum. 
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